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TECHNICAL NOTE 3734 


TURBUIENT-HEAT-TRANSFER MEASUREMENTS 
AT A MACH NUMBER OF 3*90 
■ By Maurice J. Brevoort 


SUMMARY 


Turbulent-heat-transfer measurements were obt aine d through, the use 
of an axially symmetric annular nozzle which consists of an inner, shaped 
center body and an outer cylindrical sleeve. Measurements taken along 
the outer sleeve gave essentially flat-plate results that are free from 
wall interference and comer effects for a Mach number of 3.90 and for a 

Reynolds number range of 6.3 x 10^ to 7.0 x 10 ^. The heat-transfer- 
coefficient results are slightly higher than theoretical results for a 
Mach number of 4.0 and for a ratio of inside- surface temperature to free- 
stream temperature of 4.2. The temperature-recovery factors range from 

approximately 0.89 at a Reynolds number of 6.3 x 10^ to approximately 0.86 
at a Reynolds number of 7.0 X 10^. 

INTRODUCTION 


The design of high-speed aircraft and missiles requires engine er ing 
information about heat-transfer coefficients and temperature-recovery 
factors for high speeds that extend over a wide range of Reynolds number. 
In references 1, 2, 3, and 4, local turbulent -heat -transfer measurements 
are presented for Mach numbers of 3.03, 2.06, 1.62, and 0 . 87 , respectively 

The purpose of this investigation is to extend the work of references 
1 to 4 to a Mach number of 3*90. The same type of apparatus and method 
of reducing the data was used in this investigation as employed in ref- 
erences 1 to 4. The range of Reynolds number for which measurements were 

obtained is from 6.3 X 10^ to 7.0 X 10^. The ratio of the inside-surface 
temperature to free- stream temperature varied from 4.2 to 3 . 5 . 
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SYMBOIS 

specific heat of sleeve material, Btu/(lb)(°R) 
specific heat of air at constant pressure, Btu/(lb)( 0 R) 

acceleration due to gravity, ft/sec^ 
heat-transfer coefficient, Btu/(sec)(sq ft)(°R) 
heat conductivity, Btu/(sec)(ft)(°R) 

Mach number 
Russelt number, hx/k 
Prandtl number, pcpg/k 
Reynolds nuiriber, pVx/p 

Stanton number, - = —— — 

R Pr pVCpg 

average wall temperature, °R 

effective stream air temperature at wall; same temperature 
whi ch gives a thermal potential which is independent of 
heat-transfer coefficient, °R 

stagnation temperature, °R 

inside-surface temperature of nozzle sleeve, °R 
free -stream temperature, 9R 
time, sec 

free-stream velocity, ft/sec 

specific weight of sleeve material, lb/sq. ft 

longitudinal distance along sleeve, ft (unless indicated 
otherwise) 


x 
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recovery factor. 


T e - T m 

*t - Too 


H dynamic viscosity coefficient, lb-sec/sq ft 

p free-stream density of air, slugs/ cu ft 


APPARATUS AND METHOD 


The apparatus consisted of an axially symmetric annular nozzle which 
was directly connected to the settling chamber of one of the cold-air 
blowdown jets of the Langley gas dynamics laboratory. The nozzle had a 
shaped wooden center body and two outer sleeves. One sleeve was con- 
structed of 8-inch-diameter, extra heavy, seamless, carbon-steel pipe, 
and the other sleeve was constructed of l/l6-inch stainless steel which 
was rolled into a cylinder and welded. Both of the cylinders were surface 
machined inside and outside to wall thicknesses of 0.388 inch and 
0.060 inch, respectively. The coordinates of the center body are given 
in table I. Another nozzle of similar construction, having a sleeve wall 
thickness of 0-750 inch, an Inside diameter of 11.000 inches, and an over- 
all length of 8i inches, was used for part of the test's. 

A detailed drawing of the apparatus is shown in figure 1 which gives 
the locations of the thermocouples and static-pressure orifices. Details 
of the apparatus and measuring equipment are given in references 1 and 2. 
The Mach number distribution is shown in figure 2. 

For this investigation, test runs were made for settling-chamber 
pressures of 31, 353, 395, and 4l0 Ib/sq in. gage. Excluding the first 
20 seconds, the pressures were maintained constant for each test run. 

The stagnation temperature started at essentially room temperature, 
decreased as the pip in g was cooled, and at an arbitrary time increased 
due to the addition of heat. This variation is illustrated in figure 3 
where stagnation temperature is plotted against time for a settling- 
chamber pressure of 353 Ib/sq in. gage. 

The wall temperature started approximately at the stagnation tem- 
perature and tended to approach the equilibrium temperature which was 
approximately 50° R below stagnation temperature. This variation is 
shown in figure 4 where wall temperature at station l4 is plotted 
against time for a settling-chamber pressure of 353 lb/sq in. gage. 

In figure 5 are plotted, for various times during the test run, the 
values of wall temperature against longitudinal distance along the cyl- 
inder. These values were used to deter min e the rate of change of the 
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REDUCTION OF DATA. 


The equations used in reducing the data are 


and 


T„ - T 
"r ' It - T„ 

(1) 

)> = vc 

Tw - T e 

(2) 

81 -PYV 

(3) 

Nu = St R Pr 

(» 


The method consists of selecting a recovery factor and then obtaining 
T e from equation (l). For each recovery factor that is selected, the 

corresponding quantity T„. - T e is determined and then plotted against 


the quantity 


wc 


d-*av 

dt 


(the heat input). 


The curve connecting these 


points is a straight line (eq. (2)). The true values of T_ and n 

6 *37 

are obtained when the line goes through zero. The slope of this line 
is the value of h. Figure 6 shows the values used in determining the 
recovery factor and h at station l4 for a settling-chamber pressure of 
353 lb/sq in. gage. 


The Stanton number is calculated from equation (3), and the Nusselt 
number is calculated from equation (4). The values of Frandtl number 
( 0 . 71 ) and viscosity of air were .taken from reference 5 and were based 
upon T^. The value of T„. used was that value measured 80 seconds 

after starting. The values of specific heat and specific weight of the 
sleeve material were also taken from reference 5* 
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RESULTS AND DISCUSSION 


Figure 7 shows the variation of local Nusselt number with local 
Reynolds number. The value of x used in evaluating these numbers has 
"been adjusted for x = 0 location (the effective beginning of the tur- 
bulent boundary layer) by the method of references 1 and 2. The x = 0 
locations are 4.0 inches downstream of the minimum station for a settling- 
chamber pressure of 31 lb/sq in. gage and 10.0 inches upstream of the 
minimum station for the other pressures. The Nusselt numbers were found 

to vary from 400 to 23,000 for the Reynolds number range of 6.3 X 10 5 
to 7*0 X 10?. For comparison, the curves for M = 4.0 and t H v /l 00 = 4.2 
of the Van Driest analysis (ref. 6) are shown. The average value of 
Tw/Tjo for the test results was 4.2 to 3. 5* Tbe' data were computed by 
using free-stream temperature to dete rmine the density and velocity. 

The wall temperature was used to determine the viscosity and Prandtl 
number. The data are slightly above the M = 4.0 curve. 


Figure 8 shows the variation of local temperature-recovery factor 
with local Reynolds number. The variation is from approximately 0.89 

at 6.3 X 10^ to approximately 0.86 at 'f.O x 10^. Also included for 


comparison are the curves for the recovery factor equal to 


1*1/3 


and Pr 1 / 2 . 
number. 


The wall temperature was used to determine the Prandtl 


CONCLUDING- REMARKS 


Turbulent-heat-transfer measurements that gave essentially flat- 
plate results were obtained for a Mach number of 3-90 and for a Reynolds 

number range of 6.3 X 10^ to 7*0 X 10^. The Nusselt numbers are slightly 
higher than the theoretically derived Nusselt numbers for a Mach number 
of 4.0. The temperature-recovery factors show a decrease with Reynolds 

number from approximately 0.89 at 6.3 X 10^ to approximately 0.86 at 
7-0 x 10? . 


Langley Aeronautical Laboratory, 

. National Advisory Committee for Aeronautics, 
Langley Field, Va., April 25 , 1956. 
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TABLE I.- CENTER-BODY COORDINATES 


x, in. 

Radius, in. 

x, in. 

Radius, in. 

-4.7 

2.000 

3.50 

3.0295 

-4.5 

2.020 

3-75 

2.9555 

-4.0 

2.120 

4iOO 

2.8834 

-3.0 

2.540 

4.50 

2.7473 

-2.5 ■ 

2.805 

5.00 

2.6221 

-2.0 

3-095 

5-50 

2.5067 

-1.5 

3.360 

6.00 

2.4002 

-1.0 

3.568 

7-00 

2.2093 

-.8 

3.635 

8.00 

2.o44o 

-.6 

3.690 

9.00 

1.8995 

-.4 

3.732 

10.00 

1.7734 

-.2 

3.760 

11.00 

1.6647 

0 

3.7707 

12.00 

1.5722 

-25 

3-7640 

13.00 

1.4950 

•50 

3.7430 

14,00 

1.4324 

•75 

3.7116 

15.00 

1.3838 

1.00 

3.6739 

16.00 

1.3486 

1.25 

3.6289 

17.00 

1.3258 

1.50 

3.5774 

18.00 

1.3139 

1.75 

3.5203 

18.782 

1.3U2 

2.00 

3.4585 

19.00 

1.3029 

2.25 

3.3927 

20.00 

I.2919 

2.50 

3.3237 

25.00 

1.2369 

2.75 

3-2520 

30.00 

1.1819 

3.00 

3.1786 

35-00 

1.1269 

3.25 

3-1042 

37.50 

40.50 

1.0994 

1.9070 
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Figure 1 


Four thermocouples at stations 
0, 6, IS, 18, 24, and 30 spaced 
■ 90* apart. Ons only at eaoh 
of the other stations. Total 
number of Ihermooouplee, 38. 

One sleeve 0.388-lnch mil thlcknsss 

/ One sleeve 0.060- Inch wall thickness 





-Mahogany 
oenter body 

— 47* 


“iCur pressure orltloes at stottons 
0,4, II, 25, and 31 spaoed 9<T 
apart. One only at each of the 
yinsr stations, Total number of 
pressure orifices, 33. 


Test setup 
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Figure Variation of wall temperature with time at 
settling-chamber pressure of 353 lb/sq. in. 


station l4 for a 
gage. 
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Figure 7 .- Variation of local Nusselt number with, local. Reynolds number 
for corrected location of x = 0. Viscosity and Prandtl number deter- 
mined for wall temperature. 
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Figure 8.- Variation of local recovery factor with local Reynolds number 
for corrected location of x « 0. Viscosity determined for wall 


temperature « 





